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Abstract 

Second order perturbative corrections to electron wavefunction are calculated here 
at generalized temperature, for the first time. This calculation is important to prove 
the renormalizeability of QED through order by order cancellation of singularities at 
higher order. This renormalized wavefunction could be used to calculate the particle 
processes in the extremely hot systems such as the very early universe and the stellar 
cores. We have to re-write the second order thermal correction to electron mass in a 
convenient way to be able to calculate the wavefunction renormalization constant. A 
procedure for integrations of hot loop momenta before the cold loop momenta integra- 
tion is maintained throughout to be able to remove hot singularities in an appropriate 
way. Our results, not only includes the intermediate temperatures T ~ m (where m 
is the electron mass), the limits of high temperature T » m and low temperature 
T « m are also retrievable. A comparison is also done with the existing results. 
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1 Introduction 

Finite temperature effects are important in extremely high temperature en- 
vironments, such as in the early Universe a few seconds after the Big Bang, 
astrophysical environments etc., where they are significant enough and can not 
be ignored in comparison with the vacuum contribution. The high temperature 
and density effects in ultra-relativistic plasma need to be incorporated, for ex- 
ample, in Quantum Electrodynamics (QED) plasma, quark gluon plasma and 
the core of dense stars. More recently renewed interest in hot and dense QED 
plasmas has been generated due to the possibility of creating ultra-relativistic 
electron positron plasmas with high-intensity lasers (« lO^^W/cm^) [1-3]. Two 
opposite laser pulses hitting a thin gold foil can heat up the electrons in the foil 
up to several MeV (~ 10^° K ). 

The particles propagating in vacuum can be assumed to be the ones with 
interactions switched off. When these particles propagate through a medium, 
several kinds of interaction processes take place. This makes the properties of 
the system different from that in which all the particles are assumed to be com- 
pletely independent of each other, behaving as freely propagating bare particles. 
When dealing with finite temperature environments in QED, where the parti- 
cles propagate in statistical background at energies around the thresholds for 
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particle antiparticlc pair production, the temperature effects need to be appro- 
priately taken into account. These effects arise due to continuous electron and 
photon exchanges between particles during the physical interactions that take 
place in a heat bath containing hot particles and antiparticles. The net statisti- 
cal effects of the background electrons and photons enter the theory through the 
fermion and boson distributions respectively. Finite temperature calculations 
also provide a guideline to estimate the density corrections, through chemical 
potential effects of the background plasma, for higher order loop corrections. 

The thermal background effects are included through the radiative correc- 
tions [4,5]. Self energies and the wavefunctions of the propagating particles 
acquire temperature corrections in this environment due to exchanges of energy 
and momentum with real particles. The exact state of all these background 
particles is unknown since they continually fluctiiate between different config- 
urations. Temperatures of interest in such a situation are in the range of a 
few MeV. Thermal propagators in real time formalism comprise of tempera- 
ture dependent terms added to the particle propagators in vacuum theory [6]. 
In finite temperature electrodynamics, electric fields are screened due to such 
interactions. 

Wc use the real time formulation [7] for calculations of wavefunction rcnor- 
malization as a second order perturbative correction in a due to the ease of 
obtaining the temperature corrections as additive terms to the usual contribu- 
tion in vacuum. Here we prefer that the loops with temperature dependent 
momenta are integrated before temperature independent variables in the rele- 
vant order o? loops in QED, and therefore, review the electron self energy in 
Ref. [8,9] earlier. This makes the calculations of integrations over loop mo- 
menta much more simpler and easier to handle. The results are obtained in a 
generalized form such that intermediate temperatures T ~ m are also included 
while the ranges of high temperature T » m and low temperature T « m, 
are retrieved from them as the limiting cases. 

In literature, the ways to compute finite temperature effects on phase-space, 
vertex, mass corrections and photon emission or absorption [10-19 ] are exten- 
sively discussed. The finite temperature wave function renormalization has been 
dealt with several approaches [12-21], specifically in the context of weak decay 
rates during primordial nucleosynthesis. They agree on using finite temperature 
Dirac spinors to obtain the corresponding effective projection operator. Differ- 
ences in the spinors presented in Refs. [20] and [21] were also pointed out [22]. 
However, their results for the case of /3-decay and related processes agreed with 
the ones which can be obtained using the approaches that had already existed 
for wave function renormalization, except for the case of a scalar boson decay 
in fermion-antifermion pair. We work out the two loop corrections to the fi- 
nite temperature wave function renormalization in the generalized temperature 
framework, for the first time. Section 2 is based on the re-examined and sim- 
plified calculations of loop correction upto two orders in a that contribute to 
electron self energy in this background. The expression for the relative correc- 
tion in electron mass is redone and the wave function renormalization constant 
is calculated in section 3. Section 4 gives discussion of the results. 



2 



2 Loop corrections to electron self energy 



At the one loop level, Feynman diagrams are calculated in the usual way by 

substituting the finite temperature electron, positron and photon propagators 
in place of those in vacuum. In real time formalism, the finite temperature 
terms remain separate at order a since the terms depending on temperature 
(hot) arc additive to temperature independent (cold) terms in tlie propagators. 
Therefore, at the one loop level the hot and cold loop momenta arc integrated 
separately. 

Due to interactions with the backgroimd, the electron and positron masses 
are known to get enhanced at one- loop and higher loop levels [4-9]. The photons 
also acquire dynamically generated mass due to plasma screening effect [23, 24, 
27]. The presence of effective mass implies the fact that the propagating particles 
constantly interact with the background. The radiatively generated thermal mass 
creates a mass shift and in physical quantities this acts as a kinematical cut-off, 
e.g., while determining the production rate of particles in the heat bath. 

Higher order loop corrections are required to study the perturbative behavior 
at finite temperature. The two loop integrals comprise a combination of cold 
and hot momenta which appear due to an overlap of temperature dependent 
and temperature independent terms from the particle propagators. The loop 
integrations involve an overlap of finite and divergent terms due to which these 
become analytically more complicated. In such situations, a prefered approach 
needs to be adopted for integrating overlapping hot and cold loop momenta, 
i.e., to integrate over hot loop momenta before integrating over cold ones, even 
at the two loop level [24]. This not only helps to simplify the loop integrations 
but allows one to handle the statistical effects more appropriately. 

The problem of renormalization in finite temperature field theories is some- 
what different from that at zero temperature due to the additional hot infrared 
divergences at finite temperatures. The temperature, however, acts as a regular- 
ization parameter for the hot ultraviolet divergences. The infrared divergences 
introduced in this framework are also appropriately removable in particle de- 
cay processes via bremstrahlung emission and absorption effects [5, 25]. This 
was studied in detail at order a for all the possible ranges of temperature valid 
in QED including T m [25, 27]. The renormalization of QED was also es- 
tablished at the one loop level for all the relevant ranges in temperatures and 
chemical potential [23, 25-27] . 



Fig. 1. Two loop electron self energy diagrams and counter terms 
The electron self energy is once again calculated here at the two loop level, from 
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the point of view of renormalization, by integrating over loop momenta in the 
electron self energy diagrams in Fig. 1(a) and 1(b). The relevant counter terms 
required to cancel the divergences in these two loop self energies are included in 
Fig. 1(c) and 1(d). Removal of overlapping divergences has been even checked 
in QED up to the two loop level for T ~ m, T << to and T >> to [9]. The 
integrations over the temperature dependent momenta are re-examined here 
and wherever needed are re-done before temperature independent momentum 
variables in loops for all the ranges of temperature that are relevant in QED. 
The overlapping loops in Fig. 1(a) gives nonzero real contribution to be 
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where 



f±{n,r) = { , ] ± , ^ A ; f±is,r) = { ^ ± , ^ A , 
{n + r) [n — r) [s + r) {s — r) 

/.,(p,,) = (/±^), 

Ei± = Ei[-m(3{r + s)] ± Ei[-m/3(r - s), 

f°° dk 
Ia = 8-it / —nsik), 
Jo k 

~ -87r6(m/3), 

— = — In , 

E 3E^v 1 + v 

I.p 27r3r2 ,1-v „ , 

^ ~ 47r[^ In ^-^{ma(TO/3) - Tc(m/3)} - 36(m/3)], 
E pE 1 — V 

^ - ^[{^' - ^"i2}6(m/3) + 4r{i In + 2}{ma(m/3) - Tc(m/3)}], 

p-^ V E'^ 6 V 1 — V 

with ^; = (p^ = £), 

o(m/3) = ln(l + e"™'^), 
6(m^) = ^(-irEi(-nm/3), 

n=l 

^ p—nmlB 

n = l 

and Ei(— a;) is the error integral given by 

dt 



Ei{-x) =- r 

J X 



t 



The non vanishing real terms from loop within loop correction in Fig. 1(b) are 
also re-examined and recalculated, wherever needed, by retaining the specific 
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order of integrating loop momenta, i.e., on the integration over the variables in 
hot momenta before the cold momenta, we get 

^p(P) = "'{g^^^CiZ-H + mT J2 (-ir+''[{/3(r + s)Ei[-m/3(r + s)] 

n.r.s—1 

m (n — r) m 2 m I 

1 7.p 1 2 £7° 

(n — s) IpI 2(n — r) m 

+r/3[/i(p, 7) + m^] Ei[-m/3(r + s)]}] 

„ 00 -m/3(r-7i-s) 

"is E (-ir'e-"'"=|.fe.){|i+(-i)-i^:-^) 

+ (2 - - ") Ei[-m/3(r - n)\ 

_^-m0{r-n) ^ ^^^^ + (-1)^+'] Ei[-m/3(r - n - s)]}]}. (2) 

In Eqs. (1) and (2), the preferred order of integration not only sufficiently cases 
the calculations but the results are also simpler as compared to those in Ref. 
[9] where this preference was not realized. The electron self energies at the two 
loop level in Fig 1(a) and 1(b) are then combined and rciarranged to obtain the 
temperature corrections to the electron mass and wavefunction. 



3 The Wavefunction Renormalization 

To incorporate finite temperature effecst on physical processes beyond the tree 
level, one needs to have a consistent method of temperature dependent renor- 
malization. As already mentioned, renormalizability of the electron mass was 
done through the order by order cancellation of singularities up to two loop level. 
It can be easily checked that the second order in a correction is much smaller 
than the first order contribution so that the perturbative behavior is valid. In 
a background, with T 0, the Lorentz invariance is broken and momentum 
independent renormalization constant is no longer sufficient. Donoghue and 
Holstein used the temperature dependent propagator to modify electron mass 
as well as the spinors accordingly [5]. 

The shift in the electron mass due to finite temperature effects is calculated 
here from Eqs. (1) and (2). For this all the finite terms in electron self energy 
upto second order in a are put together. Following Ref. [5] the physical mass 
of the electron at one loop was obtained in Ref. [25] in generalized form, by 
writing 

E(p) = A{p)E-/^ - B{p)p.i - dp), 
where A{p), B{p), and C{p) are the relevant coefficients. Taking the inverse of 



6 



the propagator with momentum and mass term separated as 
S-\p) = (1 - A)Ej- - (1 - B)p.7 - (m - C), 

and the physical mass rUphy = m + 6m^^^ + 5m^'^\ was deduced by locating 

the pole of the propagator t'^^T^^^ ■ 5m^^'' and dm^"^^ is the shift in electron 
mass due to temperature eSects at one and two loop level respectively. Using 
the same procedure, the relative shift in electron mass at the two loop level was 
obtained [9]. This is recalculated here, wherever required, and after recombining 
the similar summations it becomes: 
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Prom the reviewed expression for the electron self energy obtained in Eqs. (1) 
and (2) the relation for the wave function renormalization constant is derived. 
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This comes out to be 
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From this expression for Z^^, not only the behavior at intermediate tempera- 
tures T m can be extracted but the ranges of high temperature T » m, low 
temperature T «m, can be also retrieved from it as limiting cases. 



4 Results and Discussion 

With the preferred order of integration of hot loops before the cold ones, the 
previously calculated self-mass correction terms are redone, wherever required, 
for all the possible ranges of temperature. This preference simplifies the calcu- 
lations since the statistical effects arc taken care of through hot loop momenta 
integrations here, before the zero temperature integration variables are dealt 
with. Therefore, we have re-written the electron self energy expressions in QED 
at the two loop level that were presented in Rcf. [9]. This led to the modfied 
expression for the relative change in the electron mass at the two loop level in 
Eq. (3). From these corrections one can then retrieve the results for all tem- 
perature ranges of interest here, classified as, the high temperature T » m 
(having m(3 — > with e"™^ falling off exponentially as compared to 
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the low temperature T « m (with fermions contribution negligible) and the 
intermediate temperatures T m (by taking mj3 — >■ 1). 

We calculated here, for the first time in this generalized form, the wavefunc- 
tion renormalization constant using the; thermal contributions to the second 
order self energy diagrams. The divergences get cancelled as usual by including 
the counter terms in Fig. 1(c) and 1(d), as already checked [9]. The fermions 
do not pick any contribution from the heat bath at low temperature. Therefore, 
the second order in a corrections to the electron self energy at low temperature 
can be retrieved as a limiting case that contains contribution from hot photons 
only giving: 



T<<m Ol' 



to?*''-"'' 



The wave function renormalization constant upto two loops at T << m from 
Eq. (4) is therefore 

which is the same as that in Ref. [28]. The high temperature limit for this 
constant gives: 

. T»,n ^ ^21 A W3 \ AttT"^. 

n,r,5=l ' ' 

+ hip, 7) {/_ (n, r) ^ - /_ (5, r)} + /+ (n, r) ^ ^ 

+ {^Mp,7)-^|}/-(s,r)]}]. (6) 
m |p| 

It can be seen from Eq.(7) that the leading contribution in this range of 
temperature, T » m, at the two loop level is The two loop fermion self 

energy in QED has been calculated in detail recently [29] using the hard ther- 
mal loop resummation introduced by Braaten and Pisarski [30]. As far as the 
renormalization is concerned, the rcsummed hard thermal loops (HTL) do not 
affect it [29]. Hence Z2 does not get any contribution from HTL here. Moreover, 
it is worth noticing that the thermal corrections, at second order in a, to the 
wavefunction renormalization constant at extreme temperatures { T << m and 
T » m ) are still proportional to ^ as in case of the selfmass of electron. 
However, the expression for the intermediate temperatures is significantly differ- 
ent from the self-energy expression obtained earlier [9] . The calculations around 
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T ^ m for the mass and wavcfunction rcnormalizations arc still cumbersome 
at the two-loop level, but much less difficult than those in Ref. [9]. The renor- 
malizability of the theory at finite temperature can be explicitly checked and 
holds through the order by order cancellation of singularities. This provides a 
platform to include the general effects due to chemical potential in the hot and 
dense background, later on. With the experience of including chemical potential 
at one loop level, in real time formulation [23,26,27], it is foreseen that two loop 
self energies will be much more complicated but is still worth-doing to develop 
a calculational technique for high density hot plasmas or even superfluids inside 
the cores of neutron stars [31]. The modified wavcfunction is expected to modify 
the finite temperature contributions to electroweak processes [32] as well as the 
neutrino magnetic moments [33] up to the two loop level. 
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Figure caption 
Fig. 1 Two loop electron self energy diagrams and counter terms 
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